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Abstract 

We study the perturbative integrability of the planar sector of a massive SU(iV) matrix 
quantum mechanical theory with global SO (6) invariance and Yang-Mills-like interaction. 
This model arises as a consistent truncation of maximally supersymmetric Yang-Mills the- 
ory on a three-sphere to the lowest modes of the scalar fields. In fact, our studies mimic 
the current investigations concerning the integrability properties of this gauge theory. Like 
in the field theory we can prove the planar integrability of the SO (6) model at first pertur- 
bative order. At higher orders we restrict ourselves to the widely studied SU(2) subsector 
spanned by two complexified scalar fields of the theory. We show that our toy model satis- 
fies all commonly studied integrability requirements such as degeneracies in the spectrum, 
existence of conserved charges and factorized scattering up to third perturbative order. 
These are the same qualitative features as the ones found in super Yang-Mills theory, 
which were enough to conjecture the all-loop integrability of that theory. For the SO (6) 
model, however, we show that these properties are not sufficient to predict higher loop in- 
tegrability. In fact, we explicitly demonstrate the breakdown of perturbative integrability 
at fourth order. 
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1 Introduction 

The existence of integrable structures in large N gauge theories and superstring theory, that 
were found and investigated intensively in recent years, have led to an enormous progress in 
the understanding and the verification of the AdS/CFT correspondence pQ 1 . Above all there 
are the Bethe ansatz techniques which represent novel methods for computing the spectra on 
either side of the correspondence, i.e. the energy spectrum of string states on the AdS-side 
and the spectrum of conformal dimensions on the CFT-side. Hence, integrability allows for 
new tests of the AdS/CFT correspondence by making new data available. Moreover, one can 
even directly compare the integrable structures and the algebraic equations which encode the 
spectrum, i.e. an explicit computation of the spectrum can actually be evaded. The performed 
comparisons have predominantly confirmed the AdS/CFT conjecture, although some discrep- 
ancies still need to be resolved. 

In this article we concentrate on the integrable properties discovered on the gauge theory side. 

Let us therefore briefly recall the findings concerning J\f = 4 superconformal SU(iV) Yang-Mills 

theory (SYM), which are comprehensively reviewed in The statement is that the spectrum 

of anomalous dimensions of SYM in the 't Hooft large iV limit equals the energy spectrum of an 

integrable spin-chain system [1J[S]. According to this correspondence, single-trace operators of 

SYM are considered as spin-chain states, cf. Fig. ^ and the planar part of the SYM dilatation 

1 For details about integrability in the AdS/CFT correspondence we would like to refer to the reviews I2JGI 
and the references therein. 
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Figure 1: Spin-chain and magnon picture. Single-trace SYM operators are in one-to-one 
correspondence with translationally invariant spin-chain states. Different SYM fields corre- 
spond to different spin-alignments or different pseudo-particles called magnons. The highest 
weight field (Z) corresponds to spin up or an empty site. 

operator D becomes the spin-chain Hamiltonian 



Herein A = 5 Y m-^ is the 't Hooft coupling constant of SYM and Do is the planar dilatation 
operator of the free theory which measures the bare conformal dimensions. The characteristic 
feature of the integrability of the spin-chain system is the existence of an infinite tower of 
higher charges Q r >3 which commute with each other and with the Hamiltonian 2 : 



The dilatation operator is not known exactly but must be determined in perturbation 
theory in the coupling constant A jS]. It can be extracted from the logarithmically divergent 
part of the two-point functions which themselves are computed perturbatively. The application 
of the dilatation operator to a single-trace SYM operator, alias spin-chain state, corresponds 
to the attachment of some effective vertex. In the planar limit, this effective vertex is to be 
connected in all possible ways to a certain number of adjacent fields, cf. Fig. [21 The number 
of involved fields is called the range of the planar dilatation operator. It grows with the 
perturbative order. At first order in A the dilatation operator has range two, i.e. the associated 
spin-chain Hamiltonian consists only of nearest neighbor interactions. The higher perturbative 
contributions to the dilatation operator then cause long-range deformations of the spin-chain 
Hamiltonian. Notice that due to the extra factor of j in the loop counting in SYM and 

the spin-chain picture is shifted by one unit, in particular the one- loop dilatation operator 
corresponds to the zeroth order (i.e. A = 0) spin-chain Hamiltonian. 

Likewise to the dilatation operator, all higher charges also receive perturbative correction 

2 We do not consider the first charge Qi here. It is usually defined as the spin-chain momentum operator 
which generates a shift of all spins by one lattice site. However, as the spin-chain states originate from single- 
trace SYM operators which are invariant under cyclic permutations of the constituent fields, the momentum 
operator Qi acts as a c-number on all considered spin-chain states. 




(1.1) 



[Qr,Qs] = o 



Vr, s = 2, 3, . . . . 



(1.2) 
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Figure 2: Action of dilatation operator. This graphics shows the application of a sample 
dilatation operator of range two to some particular |in)-state. The result is a sum of |out)- 
states which are obtained by connecting the dilatation operator to all pairs of adjacent fields. 
In this example the dilatation operator is a mere permutation, and every local application 
yields just one contribution to the out-state. Generically, the dilatation operator will produce 
a sum of terms from each local application. 

of growing range. The general expansion is 

oo 

Qr = J2 XW2Qr > k fOT r = 2 > 3 >---, (1-3) 

fc=0 

where the range of Q T) k depends essentially linearly on r and k. Since the charges are only 
known up to some perturbative order, we also deal with integrability in a perturbative sense. 
E.g. when cut off at a certain order in A, the charges do not commute exactly with each other 
but only up to contributions of higher order. Only the lowest order is special in this respect 
since here A = and the commutation of the free charges is exact. 

Let us also recall how a rigorous proof of integrability would proceed. The reasoning is called 
algebraic Bethe ansatz and is reviewed in [J]. The main ingredient is the monodromy matrix, 
or T-matrix, which is a spin-chain operator depending on a spectral parameter u. The T- 
matrix satisfies a so-called Yang-Baxter equation, which implies that the trace tr T commutes 
for different values of the spectral parameter with itself: 

[trr(w) J trT(«')] = 0. (1.4) 

This fact qualifies trT(it) as a generating function for an infinite set of commuting charges, 
which can be obtained by a series expansion of this operator in the spectral parameter. 

If one of the commuting charges, or a combination of them, coincides with the spin-chain 
Hamiltonian, then the integrability of the corresponding system is proven. Beyond providing 
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the charges, the algebraic Bethe ansatz also yields algebraic equations, the Bethe equations, 
which determine the eigenvalues of the charges — in particular the eigenvalues of Q2, i.e. the 
energy spectrum. 

In practice, however, performing this kind of proof might be rather involved. The difficulty lies 
in finding the appropriate T-matrix and the way to extract the spin-chain Hamiltonian from 
it. There is no canonical prescription for doing that and even the Yang-Baxter equation is not 
unique. Therefore, one sometimes tries to find higher commuting charges by hand. In fact, there 
are arguments that the existence of the first higher charge Q3 already implies integrability [S] . 
And yet before searching for commuting charges, one should have a look at the spectrum. 
A good indication for the existence of these charges is the presence of degeneracies in the 
spectrum which are not due to some obvious symmetries of the system jH]. 

Alternatively, one can also utilize a coordinate Bethe ansatz, which is nicely reviewed in [S]. 
It is an ansatz for the wave function describing an eigenstate of the spin-chain Hamiltonian. 
The physical picture behind this ansatz is the propagation of the magnons (see Fig. ^) along 
the spin-chain and their scattering. The scattering is encoded in an "S-matrix", which is a 
scalar function of the momenta of the magnons involved in a scattering process. The system 
is integrable if the multi-particle S-matrix factorizes into a product of two-particle S-matrices. 
This point of view was stressed and discussed in ^Oj . 

The current status of evidence for integrability in SYM is the following. At first order in per- 
turbation theory it has been proven by means of an algebraic Bethe ansatz that the planar 
SYM dilatation operator [TI^ is an su(4|2, 2) integrable super spin-chain with nearest neighbor 
interactions [3]. At higher loop orders a similar proof is not known yet. However, the degenera- 
cies in the spectrum, a number of commuting charges and the Bethe equations have been found 
in different subsectors of the theory: In the su(3|2) subsector, the dilatation operator has been 
constructed up to third order on the basis of the symmetry algebra, some basic facts about 
Feynman diagrams and the BMN scaling behavior of the eigenvalues |12j . BMN scaling is a 
property predicted by the dual string theory and verified for SYM up to three loops in [T5] . 
In the 5u(2) subsector it was furthermore possible to map the planar three-loop dilatation 
operator to a known long-range integrable spin-chain [Tl], the Inozemtsev spin-chain |15j . But 
if this equivalence between the SYM dilatation operator and the Inozemtsev spin-chain was 
supposed to hold also at higher loops, then BMN scaling must be broken in the gauge the- 
ory ^1] . Therefore one has studied the mutual influence and the compatibility of integrability 
and BMN scaling, and Beisert, Dippel and Staudacher (BDS) succeeded in writing down an 
extension of the su(2) spin-chain Hamiltonian up to five-loop order which obeys BMN scaling 
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and commutes with higher charges Moreover, BDS proposed asymptotic Bethe equations 
and eigenvalue formulas for all higher charges to all-loop order 3 . Up to five loops these equa- 
tions reproduce the spectrum of the explicit Hamiltonian, and beyond five loops they define 
a novel spin-chain system. As the system is specified in terms of Bethe equations, the BDS 
spin-chain is in a sense integrable by definition. However, a proof that the charges and the 
Bethe equations can be derived from a T-matrix is still missing, let alone the proof that the 
BDS spin-chain indeed describes planar SYM in the su(2) subsector. 

In a parallel development, the planar integr ability of plane- wave matrix theory (PWMT) — a 
matrix model description of M-theory on a plane- wave — was found ^Hl- PWMT is very 
closely related to SYM and concerning integrability it serves as a toy model for the field 
theory. The energy operator of PWMT, which will be defined as a similarity transformation of 
the PWMT Hamiltonian below, corresponds to the dilatation operator of SYM. As a matter of 
fact, the planar energy operator coincides with the planar dilatation operator up to three-loop 
level in the su(3|2) subsector after appropriate identification of the parameters |12M18| . Hence 
all results concerning integrability that were obtained in SYM up to third order immediately 
carry over to PWMT. But as PWMT is a quantum mechanical theory, explicit calculations 
could be pushed to fourth perturbative order in the su(2) subsector of that theory [2U]. These 
computations showed that all integrability criteria such as a degenerate spectrum, existence of 
conserved charges and factorized scattering persist within PWMT to four-loop order, whereas 
BMN scaling gets violated. It also turned out that the spin-chain associated to planar PWMT 
is neither equivalent to the BDS-spin-chain nor to the Inozemtsev spin-chain. Hence, if PWMT 
should be exactly integrable in the planar sector, it would define a further long-range inte- 
grable spin-chain system. 

In summary, there are two large N matrix theories, planar SYM and planar PWMT, which 

display integrable features to relatively high perturbative order. Moreover, there is the beautiful 

all-loop conjecture of BDS. But despite the tempting simplicity of their formulas, it is strictly 

speaking absolutely unclear why the integrability should extend to all orders or even exist 

non-perturbatively. In fact, there is a number of important open questions: What are the 

essential properties a matrix theory must possess for being integrable at large N? What is the 

mechanism for the symmetry enhancement in the planar limit? What is the simplest matrix 

model leading to a long-range integrable spin-chain? 

3 "Asymptotic" refers to the fact that these equations and formulas are valid only for states longer than a 
certain threshold which increases with the considered order in the coupling constant. The finite dimensional 
subsector of shorter states is currently still inaccessible because of the so-called wrapping problem |5lll6|. which 
got its name from the fact that it is due to interactions that wrap entirely around the states. 
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In order to address these questions we study the planar integrability of an SO (6) matrix 
model. We introduce this model in the next chapter and connect it to PWMT and SYM. Then 
we investigate the model with the same tools as used in the study of those theories. We compute 
the planar energy operator in the su(2) subsector up to fourth order. Very interestingly, we find 
that this simple toy model indeed passes all commonly required integrability checks (degenerate 
spectrum, commuting charges, factorized scattering) up to and including third order — but not 
beyond. Without premonition the degeneracies in the spectrum are lifted by the fourth order 
piece of the energy operator and therefore perturbative integrability abruptly vanishes as well. 
These findings certainly sound a note of caution also for PWMT and SYM. This example 
shows that exact integrability cannot be taken for granted even if perturbative integrability 
reaches to high loop orders. 

In the following we present our results in detail. We also review the essentials of the applied 
methods. A full account for the technical details can be found e.g. in |2U1 I21| . 



2 The SO(6) matrix model 

The degrees of freedom of the model under consideration are comprised in a time-dependent 
SO(6) vector (a = 1, . . . , 6) that takes values in the adjoint representation (m = 1, . . . , iV 2 — 1) 
of su(N): 

X?(t) . (2.1) 

It is convenient to introduce the basis elements {T m ) r s of su(iV) which are traceless, hermitian 
matrices carrying a fundamental index s and an anti-fundamental index r with s,r = 1, . . . , N. 
They satisfy the relations 

[T m ,T n ] = if mn PT p , tr T m T n = 5 mn and (T m )^(T m )^ = 5 r J l s - ^5 r X . (2.2) 

We use these SU(iV) generators to define the matrix model fields as 

X a = X™T m . (2.3) 

Now, we wish to study the model given by the action 

S = Jdt tr[\D t X a D t X a - \ (f ) 2 tr X a X a + \ tr[X a , X b ][X a , X b ]\ , (2.4) 

where Dt = dt — i[ui, ] is the covariant derivative containing the scalar gauge field uj = uj m T m , 
and M is a real parameter which sets both the mass scale of the excitations and the inverse 
interaction strength as we will discuss below. As a matter of fact, this model is a consistent 
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truncation of PWMT to its SO(6) sector By definition, this means that the equations of 
motion of PWMT are satisfied by the solutions of the equations of motion of the SO(6) model 
when setting all other PWMT fields therein to zero. As furthermore the PWMT itself is a 
consistent truncation of SYM we have that the SO(6) model is also connected to the field 
theory by such a procedure. 

We want to stress that the consistency of a truncation only implies the equivalence of the 
dynamics of the selected fields at the classical level. At the quantum level, however, the fields 
which are omitted in the reduced theory will generically contribute to the dynamics of the 
mother theory. Hence, the SO(6) model should not be confused with the SO(6) subsector of 
neither PWMT nor SYM. Even within common subsectors, these are really three different 
quantum theories, which have a connection only at the classical level 4 . The relations between 
the models and their degrees of freedom are depicted in Fig. |21 

We now fix the gauge u = 0, change to the Hamiltonian formulation, and quantize the 
model canonically. The Hamiltonian reads 

H = H + V (2.5) 

with 

H = \ tvP a P a + \ (f ) 2 trA a A a = M- ti a\a a , (2.6) 

V = -\tr[X a ,X b }[X a ,X b }. (2.7) 

Here we have already inserted the following mode expansion 

Pa = \f^(a a + a\) , X a = ^(a a -ai) (2.8) 

into the free part of the Hamiltonian Hq. The modes obey the usual oscillator algebra 

K,al n ]=5 ab 5 mn , K,a%} = [ai m ,al n }=0. (2.9) 

As a matter of course, the ground state |fj) is defined as the unique state in the kernel of all 

annihilation operators a™. The excited states, which are SU(iV) invariant and relevant in the 

large N limit, are obtained from this ground state by applying a single-trace of an arbitrary 

4 It is true that, up to third perturbative order, the planar spectra of PWMT and SYM coincide in the 
largest common and closed subsector, su(3|2), after appropriately identifying the coupling constants, but this is 
a highly non-trivial fact, which still needs to be derived from first principles. At any rate, such an equivalence 
does not exist between the SO(6) model and PWMT or SYM. 
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Figure 3: Embeddings. The starting point is the full SYM theory with infinitely many ele- 
mentary fields in the singleton representation of su(4|2, 2). Seventeen fields can be identified 
as the degrees of freedom of PWMT. The largest common subsector of SYM and PWMT which 
is quantum mechanically closed in both theories is spanned by five fields transforming under 
the fundamental representation of su(3|2). The six fields of the SO(6) matrix model represent 
another common subsector of SYM and PWMT, but one which is not closed in the quantum 
theories. The su(2) subsector spanned by two complex scalar fields, however, is closed in any 
of the three models. 

number L of the matrix creation operators a\ = a}™T m : 

|V}=tr(at i at 2 ...atj|0). (2.10) 

The number L is called the length of the state. There are no states of length one due to the 
tracelessness of a a . The free energy of (t2~TUl) is E = %L. The states have a spin-chain 

interpretation as explained in the introduction, cf. Fig. ^ 

Below, we will be particularly interested in an su(2) C so (6) subsector of the model. In 
this sector the elementary oscillators form the following su(2)-doublets 

J - - _L ( a i + ia i\ = (#\ i_ (4> l \ J_( ai -ia 2 \ f Z\ 

^ \4) V2\at + ial) AW ' * V2\a 3 -iaJ " \W J 

(2.11) 
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with commutation relations 



w m Af] =^ mn , w m ,^ n ] = [<t>rJ/ i ]=o. (2.12) 

In the magnon picture, the upper component cf>\ = represents an empty side and the lower 
component represents a magnon. In the su(2) subsector there is only one kind of 

magnon. 

3 Derivation of the associated spin-chain 

As described in the introduction, the planar limit of a matrix theory has a natural interpreta- 
tion as a spin-chain system. Technically, the spin-chain Hamiltonian is given in essence by the 
't Hooft large iV limit of the so-called energy operator. In [2U] the energy operator 

T = U~ 1 HU (3.1) 

has been defined as an operator obtained from the Hamiltonian H by a similarity transforma- 
tion with the property that it does not mix states of different free energy, i.e. 

[T,H o ]=0. (3.2) 

The primal significance of the energy operator concerns the computation of the quantum me- 
chanical corrections to the free energy spectrum. Due to (|3.1|) it possesses the same eigenvalues 
as the full Hamiltonian. However, as it has no overlap between states of different free energy — a 
consequence of (|3.2j) — the energy operator disentangles the mixing problem: only the mixing 
of states within a degenerate subspace needs to be considered, the influence from states outside 
is already taken into account in (|3.1|) . 

At this point, however, we want to study the energy operator of the SO(6) model with 
respect to planar integrability. We will use the methods developed for the corresponding in- 
vestigations in PWMT and SYM. The plan of action now is the following. At first we compute 
the energy operator T in perturbation theory. Then we define a shifted and rescaled energy 
operator 

D := ^ [T - (Q\T\0)) (3.3) 

in order to have a well-behaved 't Hooft limit. The planar part D of this operator then defines 
the spin-chain Hamiltonian Q2 similar to (jl.lj) . 

For the computation of the energy operator we use the formulas derived in |181 I21| . When 
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adopted to the model (|2.5jl . they read 

oo 

T = H , T = Y,T 2k (3.4) 

t 2 = J2 n EVn E , 



k=0 



n=Yl n E VA E vn E , 

E 

T G =Y, n E [va e va e v - VA 2 E vn E v] n E 
E 

T 8 = ^2n E [VA E VA E VA E V 



- VA E VA 2 E Vn E V - VA 2 E VA E Vn E V - VA 2 E Vn E VA E V 

+ VA%vn E vn E v]n E , 

where U E is a projector onto the subspace of states with free energy E, and A E is the "prop- 
agator" defined by 

A-=E^' (3-5) 

F^ E 

The sums are taken over all free energies which are realized in this model, i.e. E £ 4^ (N\ {1}). 
The piece is called the /c-th loop contribution. All half-loop contributions are zero in this 
model, but we maintain the way of indexing the parts of T as in PWMT and SYM for a 
better comparison. Note that (|3.2|) does not uniquely specify the energy operator. In JHj this 
ambiguity was fixed by demanding to have the least number of terms. This was an essential 
requirement to conduct the highly involved computations on current computers but came 
at the price of the non-hermiticity of the energy operator T. However, by a further similarity 
transformation on top of we will change to a hermitian energy operator later. The largest 
computational effort consists in normal ordering the expressions (j3.4[) . Up to one-loop the result 
is given by 

1 f 15(iV 3 — N) + 10iVtrat aa + i : tr [a\ , a a ][aj, a 6 ] : 



M 2 L 



ti[a\,a\)[a a ,a b ] - \ : tr[al, a b ][al, a b ] : . (3.7) 



We refrain from printing the full, non-planar higher loop contributions as they are rather 
lengthy. In the su(2) subsector, however, we will give the planar part of the energy operator 
up to fourth perturbative order, below. 
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The next step is to take the planar limit. We will briefly review how this is done at the 
operatorial level and explain the necessity for the redefinition (|3,3j) . As a basic principle, the 
't Hooft limit j^2] consists of sending the rank of the gauge group to infinity and the coupling 
constant to zero with the product of these two quantities kept fixed in such a way that physical 
quantities are neither divergent nor trivial. 

The coupling constant of the SO(6) model is given by 

° 2 ■= w < (3 - 8) 

where the factor of two has been inserted for convenience (and in analogy to PWMT). This 
can easily be seen from expressing Hq and V in terms of oscillators, see (|2.Hjl — (|2.8|) . It follows 
that 

H ~M , V~j^ , IIe~1 , A E~jj- (3-9) 

Now, either from 

VI 1 M ( 1 \ k . 

JT ~W or T2k ~W^~T{w) (3 ' 10) 

we infer that the effective coupling is ~ ■jAg. This observation justifies the perturbative treat- 
ment of the model for large M, and the 't Hooft limit will also involve M — > oo. In fact, we 
will have 

N, M -> oo with A := G 2 N = ffs = fixed . (3.11) 

Note that the parameter M does not only determine the coupling constant but at the same 
time also represents the energy scale as the global factor of 4f in Q3.1U[) shows. In order to have 
a finite limit (|3.11|) . we need to rescale the energy operator. This explains the overall factor in 
()3.3|) . Hence, D measures the energies in units of an elementary excitation. 

Beyond that, we also need to analyze the dependence on N. All factors of N originate 
from closed fundamental su(iV) index loops, tr 1 = 8^ = N. Apart from the explicit factors in 
(|3.7I) . there will arise further powers of iV when the operator is applied to a state. Hence, for 
the purpose of counting its order in N, we insert the operator between two states of sufficient 
length (to avoid wrapping, cf. footnote |31 on page EJ) and normalize with respect to the case 
where no operator is inserted. The counting is most conveniently done in double line notation. 
Tab.^shows all effective one-loop vertices in the SO(6) model sandwiched between two states. 
The most important thing to note is the fact that the non-planar graphs are suppressed by 
factors of -™ compared to the planar ones and hence can be eliminated by taking iV large. 
This outcome is nothing than a specific example of the general result of 't Hooft |22| . But 
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free contractions 
for reference 



connected 
planar non-planar 



vacuum bubble 
planar non-planar 



< > 



tr a) a) a a 
~ N 



V.' 



Ntra^a 
~ N 





< 









: tr a^aa^a 




N 3 ■ 1 

~ A^ 3 



N ■ 1 

~ AT 



Table 1: One-loop efTective vertices. If we want to count the powers of A" corresponding 
to a given operator, we insert it into the two-point function and count the additional closed 
su(N) loops as compared to two-point function without operator insertion. 



one can also observe that the graphs with vacuum bubbles possess an additional factor of 
iV 2 with respect to the connected graphs. These contributions can be isolated by computing 
(0|T|0), and they must be subtracted as in (|3.3j) in order to have a finite limit ()3.11j) . Physically 
this corresponds to measuring the energy shifts less the shift of the ground state energy. In 
fact, this is a reasonable quantity to compute as one can only measure excitations above the 
ground state. In PWMT, the ground state is protected by supersymmetry and we have exactly 
<0|T|0> = 0. 

Now, we can extract the planar part D of the redefined energy operator D and, still working 
at one-loop, we find 

oo 

D = tvala a , D = ^A fc D 2fc (3.12) 

D 2 = 10tra+a a + jj tr a\a\a b a a - jj tr a\ala a a b + tr a\a\a b a b . (3.13) 

This operator is to be applied to a single-trace state (|2.1U|) in a planar fashion, i.e. the two 
annihilation operators in the traces of length four have to act onto adjacent creation operators 
of the state. This application will produce a factor of which cancels the one in ()3.13|) . The 
planar action of D can actually be described much simpler by means of the identity operator I, 
the permutation operator P and the trace operator K. They are defined to act as 

\a\al = a\a\ , P a\a\ = a\a\ , K a\a\ = a\a\8 ab . (3.14) 
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We find 

L L 

D = ki+i = L > °2 = 2 [ n " 2 Pi > i+1 + Ki = i+1 ] • ( 3 ' 15 ) 

1=1 1=1 

The summation runs over the length of the state where D is applied to. In indices Qij mean 
that the corresponding operator acts onto the i-th and j-th oscillator. The position L + 1 is 
identified with the first position. 

In this form, the experienced reader will immediately recognize D2 as an integrable spin- 
chain Hamiltonian [23] • The integr ability hinges on the ratio of the coefficients of P and K. If 
the spins transform in the vector representation of SO(n), integrability requires this ratio to 
be — — 1). For the case at hand, where n = 6, this condition is satisfied. The integrability of 
this model is hence completely established in terms of an algebraic Bethe ansatz. The details 
can be found e.g. in the review [7] and also in the original paper [I] wherein integrability in 
SYM was discovered. 

Let us now proceed to higher loop orders. Their contributions are considered as long-range 
deformations of the one-loop piece. The full spin-chain Hamiltonian is "preliminarily" defined 
as 

Qr el = }(D-Do). (3.16) 

We have called this operator preliminary as we are going to make some minor but convenient 
redefinitions. 

From now on, we will concentrate on the su(2) subsector which is generated by the fields 
4>\ defined in (|2.11|) . This entails a huge simplification, as the trace operator K annihilates all 
states in this sector. That is basically because there is no invariant tensor 5{j in su(2). In the 
evaluation of the expressions (|3.4j) we may now discard all terms which contain annihilation 
operators outside the su(2) subsector. By virtue of being a closed subsector, there are con- 
sequently also no terms with creation operator outside this sector. The computation is most 
straightforward but very cumbersome due to the plethora of terms generated in the process of 
normal ordering. We have used a [21] program which does the job for us in roughly 416 hours 
on a ordinary 2 GHz computer. The program code is printed in Add. lAl 

In order to write the result in a compact form, we adopt the frequently used notation for 
multi-permutation operators 

L 

{ni,n 2 , . . . ,ni} : = ^ P {} := L , (3.17) 
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which was firstly introduced in |B]. There are some obvious relations for these operators 



{. . . ,n,n, ...} = {...,...} , 

{. . . , n, m, . . .} = {..., m, n, . . .} , for \n — m\ > 2 (3.18) 

{ni,n 2 , • • •} = {ni + m, n 2 + m, . . .} , 



and another one which is a specialty of the su(2) sector 

{...,...} + {. .. ,n± l,n, ...} + {...,n,n± 1, ...} 
- {... ,n, ...}-{... ,n± 1,. ..}-{... ,n,n± l,n, ...} = . 



(3.19) 



Now, we can write down the planar su(2) energy operator of the SO(6) model up to fourth 
perturbative order as 

oo 

Q 2 ,o = 2{} - 2{0} , Q 2 = A fc Q 2 ,2fe (3.20) 

Q 2 ,2 = -45{} + 49{0} - 2({0, 1} + {1, 0}) , k= ° 
Q2A = ^{} - ^{0} + H6({0, 1} + {1, 0}) 

+ 4{0,2}-4({0,l,2}+{2,l,0}) , 
Q 2>6 = ~™{} + «{0} - ^P({0, 1} + {1, 0}) - ^{0, 2} 

+ S§5 ({0, 1, 2} + {2, 1, 0}) - §({0, 2, 1} + 220{1, 0, 2}) - § {1, 0, 2, 1} 

- 4{0, 3} + 4({0, 1, 3} + {0, 3, 2}) + 4({1, 0, 3} + {0, 2, 3}) 

- 10({0, 1, 2, 3} + {3, 2, 1, 0}) + f|({0, 1, 3, 2} + {2, 1, 0, 3}) 

- ^({0, 2, 1,3} + {1, 0, 3, 2} + {1,0, 2, 3} + {0, 3, 2, 1}) . 

It follows from Q^ cl through 

Q 2 = W- X Q%*W - Q\ . (3.21) 
The similarity transformation by means of the operator 

W(A) = e AAl e A2A *e A3A3 , (3.22) 

where A x = |{0}, A 2 = § ({0, 1} + {0, 1}), and A 3 = ({0, 1} + {0, 1}), corresponds to a 
change of basis that makes the spin-chain Hamiltonian hermitian. In this notation hermitic- 
ity corresponds to the invariance under {n\,. . . ,n{\ \—* {nj,...,ni}. Furthermore we have 
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subtracted a term proportional to the length operator: 

( 615 . 39123 a2 37226069 ,,\ ri 
Q\ = (9 - — A + g A - A 3 J • {} • (3.23) 

This operator is determined such that the sum of all coefficients in Q2 vanishes separately 
at any order. As a consequence Q2 annihilates the states tr(Z^) L |o). The reason for this 
subtraction is to take away the trivial contribution to the energy proportional to the spin- 
chain length from the contribution that originates purely from the magnons W*. Of course, 
the redefinitions 1)3.21(1 do not harm integrability; whenever Q2 is an integrable spin-chain 
Hamiltonian then Q^" is one as well and vice versa. Hence we will concentrate on Q2 in the 
following. 



4 Perturbative integrability and its breakdown 

In the previous section we have seen that the planar one- loop energy operator of the SO (6) 
model taken by itself defines an exactly integrable spin-chain system with nearest neighbor 
interactions. Higher loop integrability is, however, hard to prove because it is not known how to 
generalize the monodromy matrix and perhaps also the Yang-Baxter equation appropriately, 
not to mention how to show that the Hamiltonian is among the commuting charges. This is 
precisely the same situation as currently in PWMT and SYM. In those cases one therefore 
concentrates on the symptoms of integrability. These are the degeneracies in the spectra due 
to the existence of higher charges and the factorization of the S-matrix describing the multi- 
magnon scattering. These properties are widely accepted as strong evidence for higher loop 
integrability. 

In this section we demonstrate on the one hand that the SO (6) model is perturbatively 
integrable in this sense up to and including third order. On the other hand we also show 
that at four-loop level the situation is completely changed. A charge that commutes pertur- 
batively with the Hamiltonian does not exist any more and the degeneracies are indeed lifted. 
A two-magnon S-matrix does, of course, still exist but it does not reproduce the energies of 
multi-magnon states by means of Bethe equations. 

The computation of the spectrum of the spin-chain Hamiltonian (|3.20j) proceed as follows. 
First of all we need to find which states are realized and how they are organized in multiplets. 
An su(2) multiplet can be labeled by the length L and the magnon number M of the highest 
weight state. All states of one multiplet have the same length but descendent states possess 
an increased number of magnons. The multiplicities m of a certain multiplet (L, M) in the 
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L 


4 


5 


6 


7 


8 


9 


M 


2 


2 


2 


3 


2 


3 


2 


3 


4 


2 


3 


4 


m p 


1+ 


1" 


2+ 


1" 


2~ 


1± 


3+ 


1^1" 


3+ 


3" 


3± 


1±,2" 



Table 2: States in su(2) sector. Multiplets are labeled by the length L, the magnon number 
M of the highest weight states and the parity p. This table lists the multiplicities m of all 
irreducible multiplets with L < 9. Here denotes m pairs of multiplets whose partners have 
opposite parity. These pairs, printed in bold face, have degenerate energies if the associated 
spin-chain system is integrable. 

spectrum is given by the number of linearly independent M-magnon states of length L that 
are highest weight states, i.e. that are annihilated by the raising operator J + . J + acts on states 
according to the Leibniz rule and then on each single oscillators as 

J+W ] = Zt , J+Zt = . (4.1) 

The multiplets can furthermore be disentangled such that their states possess definite par- 
ity p. The parity conjugation operator P acts on the SU(iV) generators as transposition and 
multiplication by —1, i.e. on single-traces as 

P trT rni ■■■T niL P~ l = (-l) L tvT mL ■■■T mi , (4.2) 

and the vacuum has positive parity f|0) = |o). Hence, parity conjugation essentially reverses 
the oscillators in a state 

P tr ^ • • ■ 4 1 ) = {-if tr 4 • • • 4 10> . (4.3) 

The multiplicities of all su(2) multiplets with L < 9 together with their parity are listed in 
Tab. |2 Of particular importance for integrability are the parity pairs, denoted by in 
the table. The states of a pair, |+) and | — ), are related by a charge Qs which is parity odd, 
PQ3P 1 = —Q3, and which commutes with the Hamiltonian, [Q2>Qs] = 0. This implies that 
these states have identical energy. In formulas this is 

\ + )=Qs\-) , Q 2 \±) = \±)qt => 4 = <h- (4-4) 

The systematic occurrence of these degeneracies are a strong indication of the existence 
of such a commuting charge Q3 and thus also of the presence of integrability. In Tab. 01 we 
give the energies q 2 for all parity pairs up to L < 9. The observation is that the partners of a 
pair have exactly identical energies up to third perturbative order but slightly different ones 
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L 


M 


Spin-chain energy 


Bethe momenta 


7 


2 


92 


- 4 80A + 2595A 2 123 1 1 2 451 A 3 ( ) - Bethe 






pi = f 

P2 = ~Pl 






92 


- 12 264A + 9111A 2 1502247 A 3 ( )- Bethe 






Ttt 

Pl = f 

P2 = "Pl 


7 


3 


92 


f 7205509 A 3 ( , 1 

-in 21 5A 4- 29325 A 2 4- J 24 1 j 

10 215A+ 4 A +| ( } _ Bcthc 






Pi = ±1.16±0.93i 
p 2 = ±1.16 =F 0.93i 
P3 = -pi - P2 


8 


3 


92 


I- 1367125 A 3 ( +) _ Bethe 

= 8-168A±5633A 2 ± j r J 116 ^ J j 






pi = ±0.96 ± 0.59i 
p 2 = ±0.96 =F 0.59i 
P3 = -pi - P2 


9 


3 


92 


« 6.45322 - 132.568A + 4378.58A 2 + / - 175265 ' 3A ' 

[-175267.0A 3 


(+) 
(-) 


= Bethe 


pi = T 0.83 =F 0.43i 
p 2 = T 0.83 ± 0.43i 
P3 = -pi - P2 






92 


, , f -333240. 5 A 3 
f» 11.0399 - 237.742A + 8122.41A 2 + { 

[ -333230. 2A 3 


(+) 

(-) 


= Bethe 


pi = =Fl-28 =F l-26i 
p 2 = =Fl-28±1.26i 
P3 = -pi - P2 






92 


f -511941. 8A 3 
« 16.5068 - 360.690A + 12421. 3 A 2 + I 

[-511944.4A 3 


(+) 

(-) 


= Bethe 


pi = ±2.98 
p 2 = ±1.15 

P3 = "Pl - P2 


9 


4 


92 


f-=3A 3 (+) 
-10 215A± 29405 A 2 ± { 797 2 32139 A 3 ( ) 

[ 79734083 A 3 Bethe 






Pi = ±2.63 

p 2 = T 0.53 ± 0.88i 

P3 = Tl-77 

P4 = "Pl - P2 - P3 



Table 3: Selected states of su(2) spectrum. We list the eigenvalues as computed by direct 
application of the Hamiltonian to states and explicit diagonalization. "(+/-)" refers to the 
parity of the states. We observe the degeneracy of the parity pairs up to three-loop and the 
lift of the degeneracy at four-loop order. We also calculate the eigenvalues by means of the 
Bethe ansatz. (The values pi give the A -th order of the quasi- momenta.) The result coincides 
with the elementary calculation only up to three loops. 

at fourth order. 

We show that this behavior is not just an effect for short lengths L by presenting a local, 
parity-odd charge Q3 satisfying 

[Q2,Q3] = C(A 3 ) (i.e. zero up to fourth order terms) . (4.5) 
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It reads 



2 



Q3,0 



C3,0,l({0,l} 



{1,0}) 



Q 3 = ^A fe g3,2fc 



(4.6) 



fc=0 



Q 3 ,2 = C 3 ,2,l({0, 1} - {1,0}) + 2c 3 ,o,i({0, 1,2} - {2, 1,0}) , 

Q 3 ,4 = C 3 ,4,l({0, 1} - {1,0}) + (-57C8.0.1 + 203,2,1 - 2C 3 ,4,2)({0, 1, 2} - {2, 1, 0}) 

- 2c3, ,i({0, 1,3} + {0, 2, 3} - {0, 3, 2} - {1, 0, 3}) + c 3 ,4,2({0, 1, 2, 3} - {3, 2, 1,0}) 
+ (4c 3 ,o,i - c 3 ,4, 2 )({0, 2, 1, 3} - {1, 0, 3, 2} + {0, 1,3, 2} - {0, 3, 2, 1} 

+ {1,0,2,3} -{2,1,0,3}) . 



In this notation parity conjugation acts as P{n\, . . . , n;}P _1 = {— n\, . . . , — nj}. The con- 
stants c are not fixed by solely demanding the commutation with Qi- And as to be expected 
from the lift of the degeneracies, there is no parity-odd operator Q^^o of (maximal) range six 
which satisfies 



Finally, we consider the perturbative asymptotic Bethe ansatz |10| for this model. We deduce 
the two-magnon S-matrix from the spin-chain Hamiltonian Q2 and compute from it the energies 
of all states in Tab. |5] by means of Bethe equations. This will show the factorization of the 
S-matrix up to three-loop level, but it will not reveal any suspicious behavior of the system 
that might help to anticipate the breakdown of integrability at fourth order. 

We derive the S-matrix according to the method explicitly explained in |201 I21j . One 
considers all spin-chain fragments (i.e. no su(iV) trace) of length L with 2 magnons at positions 
l\ and I2 



O hh ) := {Z^) h - l W\Z^) h - h - l W\Z^) L - h \ti) for 1 < h < l 2 < L- 1 . (4.8) 



[Q3,6, C?2,o] + [Q3,4, Q 2 , 2 ] + [Q3,2, Q2A + [Q 3 ,0, Q2,G] = . 



(4.7) 



These fragments are superposed to energy eigenstates of the Hamiltonian 




(4.9) 
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according to the perturbative asymptotic Bethe ansatz as 

L 



\puP2) ■= ^2 a(h,l 2 ,Pi,P2)\Oi u i 2 ) (4.10) 



h,h=l 



with 



a(h,h,Pi,P2) = e^+P^fih -h,Pi,P2) 

+ j<*h-H»h) f ( L -l 2 + h,pi,P2)S(p 2 ,Pi) 



(4.11) 



The exponentials describe the free propagation of the magnons along the spin-chain with 
quasi- momenta p\ and p 2 ■ The functions 

f{l,Pi,P2) = 1 + A'/ (/,Pi,p 2 ) + A l+1 fi(l, P i,P2) + A l+2 f 2 (l, Pl ,p 2 ) + • • • (4.12) 

and 

S(pi,P2) = S ( Pl ,p 2 ) + AS 1 ( Pl ,p 2 ) + A 2 S 2 ( Pl , P2 ) + A 3 S 3 ( Pl ,p 2 ) + ... . (4.13) 

describe the magnon scattering due to the long-range (/) and the nearest neighbor interactions 
(S), respectively. Joining the end of the fragments by taking the trace in order to obtain cyclic 
spin-chains leads to the two-magnon Bethe equations 

exp(iLpi) = S(pi,p 2 ) , exp(iLp 2 ) = S(p 2 ,pi) (4.14) 

and the total momentum condition 

Pl +p 2 = 2irZ. (4.15) 

In an integrable system where the multi-magnon S-matrix is factorized into products of the 
two-magnon S-matrix, these equations have straightforward generalizations to the M-magnon 
case 

M M 

exp(iL Pk ) = 11 S( Pk , Pj ) for k = 1, . . . , M , = 2ttZ . (4.16) 

j=i i=i 

Specializing to the SO (6) model, we act with the spin-chain Hamiltonian (|3.20|) onto (|4.10|) 
and demanding (|4.9|) fixes all functions q 2 (pi,p 2 ), f(l,Pi,P2), and S(pi,p 2 ) up to 0(A 3 ) (fourth 
perturbative order). The energy is given by the sum q 2 (pi,p 2 ) = Q2(pi) +92(^2) of a one-magnon 
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energy 



q 2 (p) = 8sin 2 (|) - A [l48sin 2 (|) + 32sin 4 (|)] 

+ A 2 [4601 sin 2 (|) + 1472 sin 4 (§) + 256 sin 6 (§)] (4.17) 
_ A 3 [1066301 ^2(2) + i9np9 sin 4 ( £ } + 18816 sin 6(|) + 25 60sin 8 (2)] . 

In the generalization to M magnons, every magnon contributes a corresponding portion to the 
total energy 

M 

92 = J> 2 (Pi)- ( 4 -!8) 

i=l 

The first three orders of the S-matrix can be brought into the usual form 

where here the phase function reads 

<p(p) = icot(2)[l+8Asin 2 (|)- 220A 2 sin 2 (2)-32A 2 sin 4 (2)] . (4.20) 

We refrain from printing the fourth order piece as it is not of this or any other obvious compact 
form. Also the introduction of an exponential factor into l)4.19j) (as in PWMT [2U] or on the 
string theory side [23]) did not lead to a meaningful expression. Moreover, we do not print the 
function f(l,pi,P2) as it is not relevant for the computation of the eigenvalues. 

Now, we solve the general Bethe equations (j4.16|) with (|4.19f) for the cases of Tab. The 
lowest order quasi-momenta are given in the table. Then we compute the energy eigenvalues 
using (|4.18|) and compare them with the results from the direct diagonalization of the Hamil- 
tonian. Up to third order we find exact agreement. At fourth order the Bethe ansatz approach 
must fail as, on the one hand, the degeneracies between parity pairs are lifted but, on the 
other hand, the Bethe ansatz cannot distinguish between the two partners of a parity pair — 
essentially because q 2 (p) = Q2(—p)- However, we observe that in the considered three-magnon 
cases, the Bethe solutions coincides with one state of a pair. This is no longer the case for 
the four-magnon pair where the Bethe solution differs slightly from the true energies of both 
states. For the two-magnon states, on the other side, the Bethe ansatz still works at fourth 
order. This in not surprising as the Bethe equations (|4.14|) in this sector follow from a rigorous 
derivation; no factorization had to be assumed. 
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5 Conclusions 



At the moment the approved data concerning integrability in SYM reach up to third |12| and 
in PWMT up to fourth perturbative order |20j . In this article we have posed the question 
whether these data are sufficient to justify the belief in all-loop integrability. In fact, we have 
presented a toy model where the answer is negative. The model under consideration is the 
SO (6) matrix model, which is related to both PWMT and SYM by a consistent truncation. 
This model displays the same integrability properties as its mother theories: exact integrability 
at one-loop order, and degeneracies in the spectrum, existence of conserved charges as well 
as factorized scattering up to three-loop order. But for all that, as one proceeds to the next 
order, all of these properties and thus integrability abruptly cease to exist. 

The moral of this article is the simple warning, that despite the presence of perturbative 
integrability at low orders, one cannot take higher or all-loop integrability for granted. Clearly, 
also before this example it was evident that integrability may break down at any arbitrary 
order; just take any long-range integrable spin-chain, e.g. the Inozemtsev spin-chain, and add 
a deformation that breaks integrability at some order. The point here is that we did not use 
an artificially constructed spin-chain but started from a decent looking matrix theory which 
pretended to be integrable. 

Admittedly, the SO (6) matrix model has very much less structure than SYM and also 
PWMT. On the one hand, the additional symmetries of the latter theories might be just what is 
needed for exact integrability. From the three-loop investigations concerning the su(3|2) sector 
of SYM |12| . we know that the symmetry (together with BMN scaling) fixes the dilatation 
operator strongly enough to imply planar integrability. However, the higher the perturbative 
order the weaker are the restrictions from symmetry. So, on the other hand, the additional 
symmetries may just shift the breakdown of integrability to a higher level. Therefore it is 
essential to know how much freedom the symmetries ultimately leave and what the basic 
cause for integrability really is. 

It does not lead to any principle improvement if we could push the perturbative results by 
one or two orders (unless we find the breakdown of integrability) and we eventually need to 
verify the planar integrability non-perturbatively. In the field theory, however, this is a rather 
formidable aim since the dilatation operator is not known exactly. Therefore we propose to 
intensify the study of quantum mechanical matrix theories where the full Hamiltonian is given. 
A particularly interesting theory would be PWMT but other appropriate toy models can be 
constructed as well. In possession of the full Hamiltonian, one might be able to go without 
invoking perturbation theory and one can try to determine charges that commute exactly. An 
important prerequisite for such a program, however, would be to implement the planar limit at 
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the level of the Hamiltonian. This is a highly non-trivial task but it is of central and essential 
significance for proving integr ability in large N matrix theories. 
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A Program code 

We print the computer program which was used to normal order the energy operator T as 
defined in (|3.4[) . It is written in |24| and the required interpreter can be downloaded from 
http : / / www . nikhef . nl/ ~ form / . 

Starting point is the definition of the energy operator in lines 1^ 1171 Here, In marks the 
place where the initial state with free energy Eq will be attached, and D(p,e) stands for the 
p-th power of the propagator A£ 0+e (|3.5|) if p > 1 and for the projector IlE +e if P = — 1- 
The vertex V is expressed in terms of the matrix model field X(a,m) <-> X™ in lines I19H24I and 
subsequently in terms of modes A(a,m,+1) *-> a^" 1 , A(a,m,-1) <-> a™. The mass parameter M 
is set to unit during this computation. Now, we can evaluate the propagators and projectors 
( lines l3T]to43|) and perform the normal ordering of the expression (lines H5M51|) . In line I5UI we 
turn the oscillators into commuting functions A i— > Ac thereby giving up the ordering. From 
now on, all expressions are understood to be in normal ordered form. It remains to perform 
the gauge group algebra, cf. lines I53H6U1 Since we are only interested in the large N limit 
and since the potential V would not imply any interactions between u(l) fields, we may as 
well work with the gauge group U(iV) instead of SU(iV). This reduces the computation time 
enormously. After all contractions are carried out, we discard the U(l)-part again in line I6U1 
Eventually we truncate to the su(2) subsector, i.e. keep only the oscillators Z< , Z, W' , W . 

The output of the program is the planar energy operator in the su(2) subsector spoiled by 
terms that are sub-leading in and a dominant vacuum contribution. The vacuum bubbles 
can be identified easily as they do not carry any oscillators, and we erase them by hand. The 
sub-leading terms are eliminated when we change to the language of permutation operators 
(|3.17jl as follows. We apply the computed energy operator (without vacuum contribution) to 
some long sample states and keep only the leading order in N. Then we apply an ansatz for 
the planar energy operator in terms of permutation operators to the same sample states and 
fit the coefficients of the ansatz such that the out-states for both applications match. 
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Due to a limited hard disk capacity ( writes large intermediate results to disk) the actual 
computation was performed in several parts, i.e. the program was run various times — each 
time for a different parts of the input lines IH M171 

1 symbol N , e , u , power , I ; 

2 autodeclare index m , n ; 

3 autodeclare index a=6,b=6; 

4 function V , X , A , D , In ; 

5 cf unction Tr , Ac ; 

6 vector Zd , Z , Wd , W ; 

7 dimension N; 

8 



9 


L 


[T2] 


= + 


D( 


-1 ,0) *V*In; 




1(1 


L 


[T4] 


= + 


D( 


-1 ,0) *V*D(+1 ,0) 


*V* In ; 


11 


L 


[T6] 


= + 


D( 


-1 ,0) *V*D(+1 ,0) 


*V*D (+1 ,0) *V*In 


12 








D( 


-1 ,0) *V*D(+2 ,0) 


*V*D(-1 ,0) *V*In; 


13 


L 


[T8] 


= + 


D( 


-1 ,0) *V*D(+1 ,0) 


*V*D(+l,0)*V*D(+l,0)*V*In 


14 








D( 


-1 ,0) *V*D(+1 ,0) 


*V*D(+2,0)*V*D(-l,0)*V*In 


15 








D( 


-1 ,0) *V*D(+2 ,0) 


*V*D(+l,0)*V*D(-l,0)*V*In 


16 








D( 


-1 ,0) *V*D(+2 ,0) 


*V*D(-l,0)*V*D(+l,0)*V*In 


17 






+ 


D( 


-1 ,0) *V*D (+3,0) 


*V*D(-l,0)*V*D(-l,0)*V*In; 
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19 #message Substitute vertex 

20 #do x = l , 4 

21 id.once V = l/2*( + X ( a ' x ' , ml ' x ' ) *X ( a 1 x ' , m2 ' x ' ) *X (b ' x ' , m3 ' x ' ) *X (b ' x ' , m4 1 x ' ) 

22 - X (a' x ' ,ml 'x ' ) *X(b 1 x ' , m2 ' x ' ) *X (a'x ' , m3 'x ' ) *X(b ' x ' , m4 1 x ' ) 

23 ) *Tr (ml ' x ' ,m2 ' x ' , m3 ' x ' , m4 1 x ' ) ; 

24 #enddo 

25 

26 #message Replace fields X by oscillators A 

27 id X(a?,m?) = I * A ( a , m , -1 ) - 1 * A ( a , m , + 1 ) ; 

28 Id I~2 = -1; 

29 

30 #message Remove propagators and projectors 

31 repeat ; 

32 * Commute propagators and projectors to the right 



33 id D (power ? , e?) * A (a? , m? , u?) = A (a , m , u) *D ( power , e -u/2) ; 

34 * Evaluate propagators 

35 id D(l,0)*ln = 0; 

36 id D(2,0)*In = 0; 

37 id D(3,0)*In = 0; 

38 id D(4,0)*In = 0; 

39 id D (power ?{>=0} , e? ! {> = , < = 0}) *In = In *( 1 / e ) " power ; 

40 * Evaluate projectors 

41 id D(-l,0)*In = In; 

42 id D(-l ,e?!{>=0 , <=0}) *In = 0; 

43 endrepeat ; 

44 

45 #message Normal ordering 

46 repeat ; 

47 id A (a? ,m? , -1) *A (b? , n? , 1) = A (b , n , 1) * A (a , m , -1) + d_ (a , b) *d_ (m , n) ; 

48 endrepeat ; 
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49 
50 
Bl 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
fifi 
67 



* give up the ordering , normal ordering implied 
id A (?m) = Ac (?m) ; 
id In = 1; 



#message U(N) Group algebra 
repeat ; 

id Tr(?i,m?,?j)*Tr(?k,m?,?l) 
id Tr (?i ,m? ,?j ,m? ,?k) 
id Tr() 
endrepeat ; 

* Discard U(l) part now 
id Tr(m?) = 0; 



Tr (?i , ?1 , ?k , ? j ) ; 
Tr (? j ) *Tr (?k ,?i) 
N; 



#message Truncate to su(2) 



id Ac (a? ,m? ,+1) *Ac (a? . 



,+D 



id Ac (a? ,m? , -1) *Ac (a? ,n? , -1) 
id Ac (a? ,m? ,+1) *Ac (a? ,n? , -1) 

cycl Tr; 



= 0; 
= 0; 

= Zd(m)*Z(n) + Wd(m)*W(n); 



69 print +s +f ; 

70 

71 . end 
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